The crystalline electric field as a probe for long range antiferromagnetic order and 

superconductivity in CeFeAsOi_a;Fa. 
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We use inelastic neutron scattering to study the crystalline electric field (CEF) excitations of 
Ce^+ in CeFeAsOi-a;Fi,(a; = 0,0.16). For nonsuperconducting CeFeAsO, the Ce CEF levels have 
three magnetic doublets in the paramagnetic state, but these doublets split into six singlets when 
Fe ions order antiferromagnetically. For superconducting CeFeAsO0.84F0.i6 {Tc = 41 K), where the 
static AF order is suppressed, the Ce CEF levels have three magnetic doublets at huo = 0, 18.7, 58.4 
meV at all temperatures. Careful measurements of the intrinsic linewidth F and the peak position 
of the 18.7 meV mode reveal clear anomaly at Tc, consistent with a strong enhancement of local 
magnetic susceptibility x"(^'^) below Tc- These results suggest that CEF excitations in the rare- 
earth oxypnictides can be used as a probe of spin dynamics in the nearby FeAs planes. 



The rare-earth (i?) oxypnictides with general formula 
i?FeAsO [R = La, Sm, Ce, Nd, and Pr) are currently 
attracting much attention due to the discovery of high- 
transition temperature (high-T^) superconductivity in 
these materials upon chemical doping [l|, y, y, Ij, |5|] . Al- 
though superconductivity in electron doped LaFeAsO ap- 
pears at a moderate superconducting temperature of 28 
K [l|, replacingjja with other rare earth ions increases 
Tc up to 55 K [2, S 1J|, making the rare-earth oxypnic- 
tides a new class of high-Tc superconductors with critical 
temperatures only surpassed by high-Tc copper oxides. 
Since superconductivity in these rare-earth oxypnictides 
appears after electron-doping to suppress the static an- 
tiferromagnetic (AF) order in their parent compounds 
[fil, LZI, Isl, Isl, llO| , it is important to determine the influence 
of magnetic interactions on the superconducting proper- 
ties. Compared to LaFeAsO, where Fe is the only pos- 
sible ion carrying a significant magnetic moment, the 
rare-earth oxypnictidc with unpaired 4/ electrons such 
as Ce^"*" in CeFeAsO offers an unique opportunity to 
study the interplay between the rare-earth and Fe mag- 
netic ions. In particular, by using neutron scattering to 
study the crystal-electric-field (CEF) excitations of the 
rare-earth in i?FeAsO and their doped superconductors, 
one can determine the electronic ground state of the R 
ions and therefore understand the low temperature ther- 
modynamic and magnetic properties of these materials 



Ul, ll3l- Furthermore, since the R ions are situated near 
the superconducting FcAs layer, CEF excitations at the 
R sites are sensitive to the electronic properties of the 
FeAs layer, and can be used as a probe to study the in- 
fluence of superconductivity on spin fluctuations in the 
FeAs plane [ll,[li,[ll. 

In this paper, we report inelastic neutron scatter- 
ing studies of the Ce'^"'" CEF excitations in the AF or- 



dered CeFeAsO and superconducting CeFeAsO0.84F0.i6 
(Tc = 41 K) 0. For CeFeAsO, we find that the 
Ce'^+ CEF levels are composed of three magnetic dou- 
blets at fiLo = 0, 18.7, and 60 meV in the paramag- 
netic state. When the Fe long range AF order sets in, 
these doublets are split into six singlets. In the case 
of CeFeAsO0.84F0.i6j although the three doublets are no 
longer split, the intrinsic linewidth F and the peak posi- 
tion of the TiLOi ~ 18.7 meV mode show a clear anomaly 
below Tc- These results suggest that the linewidth and 
position of the CEF transitions in the rare-earth oxyp- 
nictides are sensitive to the Fe spin ordering, and can be 
used as a direct probe of the spin dynamics in the nearby 
FeAs planes. 

Our experiments were carried out using the MER- 
LIN chopper spectrometer at ISIS facility, Rutherford- 
Applcton Laboratory, Didcot, UK [11]; the BT4 Filter 
Analyzer Neutron Spectrometer (FANS) at the NIST 
center for neutron research (NCNR) [l7[; and the 
NG4 Disk-chopper time-of-flight spectrometer (DCS) at 
NCNR, Caithersburg, Maryland. MERLIN is a high 
count rate, medium energy resolution, direct geometry 
chopper spectrometer with a large solid angle of posi- 
tion sensitive detectors. FANS is a high count instru- 
ment to measure inelastic excitations on powders and 
DCS is a cold neutron direct geometry chopper spec- 
trometer. Our CeFeAsO and CeFeAsO0.84F0.i6 sam- 
ples were prepared using methods described in Ref. [2| 
and their structural/magnetic properties are discussed in 
Ref. [71]. For MERLIN measurements, a low temperature 
He-4 cryostat and closed cycle He-gas refrigerator were 
used for the temperature variable. The incident beam 
energies were E, = 7,10,30,70,105,300 nicV. To sepa- 
rate the CEF magnetic scattering from phonon excita- 
tions, we also measured LaFeAsO and LaFeAsOo.88Fo.o8 
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FIG. 1; (Color online) Summary of the CeFeAsO crys- 
tal/magnetic structure and CEF levels determined from our 
inelastic neutron scattering experiments, (a) The Fe spin or- 
dering in the CeFeAsO chemical unit cell, (b) The Fe spins in 
CeFeAsO with respect to the Ce positions. The Fe moments 
lie in the a-b plane along the a-axis and form an antiferro- 
magnetic coUinear spin structure, (c) Ce'^"'" CEF levels in 
CeFeAsO for temperatures above and below the Fe AF Neel 
temperature of Tn ~ 140 K [J- The arrows indicate pos- 
sible transitions, (d) Ce"^""" CEF levels in superconducting 
CeFeAsO0.84F0.i6 at low temperature. 



samples [g as reference materials for phonon subtrac- 
tion from the Ce samples. The data for CeFeAsOi-a-F^; 
and LaFeAsOi_a;Fa; were converted to absolute units 
of mb/sr/mcV/f.u. by normalizing the scattering from 
vanadium measurements made with the same incident 
energies. For the DCS measurements, we used a He-4 
cryostat and an incident beam energy of 3.55 mcV and 
for the FANS measurements we used a top-loading Hc- 
gas refrigerator. 

Figure 1 (a) shows the position of the Ce ion in the crys- 
tal structure environment of CeFeAsO. Relative to the Fe 
sublattice, the Ce'^+ ions are located alternately above 
and below the (AF ordered) Fe layers as shown in Fig. 
1 (b) . Figures 1 (c) and 1 (d) summrize the Ce CEF excita- 
tion energies determined from our inelastic neutron scat- 
tering experiments for CeFeAsO and CeFeAsO0.84F0.i6j 
respectively. According to neutron powder diffraction ex- 
periments [7|, both CeFeAsO and CeFeAsO0.84F0.i6 have 
a tetragonal (space group PA/nmm) crystal structure at 
room temperature. However, on cooling, CeFeAsO first 
exhibits a structural phase transition, changing the crys- 
tal symmetry from tetragonal to orthorhombic (space 
group Cmma), and then orders antifcrromagnctically 
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FIG. 2: (Color online) Temperature dependence of the CEF 
excitations in CeFeAsO and CeFeAsO0.84F0.i6 and our model 
determination of the CEF levels, (a) Ce CEF magnetic ex- 
citations in absolute units after subtracting the LaFeAsO 
phonons. The solid line is our model fit with parameters 
listed in Table I. The inset shows the raw DCS data integrated 
from 0<(3<2.2A~^in arbitrary units as a function of tem- 
perature without the LaFeAsO phonon subtraction, (b) Ce 
CEF excitations at 200 K; solid line is our model calculation, 
(c-d) CEF excitations and their temperature dependence for 
CeFeAsO0.84F0.i6- The solid lines are model calculations. In- 
set in (c) shows that the 0.7 meV excitation seen in CeFeAsO 
is missing in CeFeAsO0.84F0.i6. Error bars represent one stan- 
dard deviation. 



with a spin structure as shown in Figs, la and lb [6J, |7[; 
CeFeAsO0.84F0.i6 maintains the tetragonal structure for 
all temperatures and does not order magnetically above 4 
K. In the tetragonal structure, the Ce atoms are located 
at the 2c crystallographic site which has C^v point sym- 
metry. This gives three non-zero CEF parameters in the 
crystal field Hamiltonian and its form in Stevens opera- 
tors formalism is HcBi^^ (041, ) = B^Ol+BlOl+BlOl. In 
the case of the low-temperature orthorhombic structure 
of CcFcAsO, the Ce atoms are at the Ag (0, 1/4, z) site 
which gives local point symmetry of mm2 {C2v)- The 
resulting Hamiltonian then involves five non-zero crys- 



tal field parameters as Hcef{C2v) = B2O2 + S|0| 



B\0\, where i?™'s are the CEF param- 



BlOl + BlOl^^,^^, 
eters to be determined from the experimental data; and 
O™ 's are operator equivalents obtained using the angular 
momentum operators 11|. 

We collected neutron scattering data on CeFeAsO 
and CeFeAsO0.84F0.i6 on MERLIN and DCS spectrome- 
ters with different incident beam energies to search for 
Ce'^+ CEF excitations. To eliminate phonon scatter- 
ing, we carried out identical scans using LaFeAsO and 
LaFeAsOo.ssFo.os as reference materials. Figures 2(a) 
and (b) show the phonon subtracted energy scans for 
CeFeAsO at 7 K and 200 K on MERLIN obtained with 
Ei = 105 meV. The inset of Fig. 2(a) plots similar data 
taken on DCS with Ei = 3.55 meV. Identical scans taken 
for CcFcAsO0.84F0.i6 are shown in Figs. 2(c) and (d). 



TABLE I: Refined B^ CEF parameters for CeFeAsO and 

CeFeAsO0.84F0.i6. 

CoFoAsOo.84F0.i6 CeFeAsO (> Tn) CeFeAsO (< Tn) 



Bl 


2.29 ±0.11 


3.19 ±0.15 


1.648 ± 0.035 


Bl 


-0.061 ±0.0064 


-0.032 ± 0.0072 


0.1073 


Bl 






-3.098 ± 0.064 


Bi 






-0.288 ±0.009 


Bi 


0.698 ± 0.023 


0.755 ± 0.049 


0.591 ±0.018 



We first discuss results on the superconducting sample as 
there are no complications of structural distortion and Fc 
magnetic order. To obtain the i3™'s CEF parameters in 
the HcEpiCiv) Hamiltonian, we first used the FOCUS 
program, which has a Monte Carlo search routine, to fit 
the observed two CEF excitations at 18.7 meV and 58.4 
meV in Figs. 2(c). We then fit many spectra simualte- 
niously with different incident energies and temperatures 
using a CEF fit program and the results are plotted as 
sofid fines in Figs. 2(c) and (d) for 7 K and 200 K, re- 
spectively. Tables I and II summarize the i3™'s CEF 
parameters and wave functions for CeFeAsO0.84F0.i6 ob- 
tained from those fits. 

Comparing to the superconducting CeFeAsO0.84F0.i6, 
the Ce CEF excitations in CcFcAsO near 19 and 65 meV 
have clear double peaks at low temperature that become 
a single peak at 200 K [Figs. 2(a) and (b)]. In addition, 
the low-energy spectra in the inset of Fig. 2(a) shows 
a clear peak around 0.7 meV that is not present at 171 
K. To understand this phenomenon, we carried out care- 
ful temperature dependence measurements of the ^19 
meV CEF excitation. Figure 3(a) shows the raw S{Q,uj) 
spectra of CeFeAsO collected on MERLIN at 60 K using 
Ei = 30 meV. After subracting the phonon scattering 
background collected using LaFeAsO, the Ce CEF level 
shows two clear bands of excitations at 16.6 and 20.4 
meV [Fig. 3(b)]. Figure 3(c) shows the detailed temper- 
ature dependence of the ~19 nieV exitations and Figure 
3(d) plots splitting of the two low temperature peaks. 
Comparison of these figures with the Ncel ordering tem- 
perature of CeFeAsO in Figure 3(e) makes it immediately 
clear that the CEF splitting is due to the long range AF 
Fe ordering [7|. 

In principle, the orthorhombic structural distortion 
that precedes the AF ordering in CeFeAsO can have an 
effect on the Ce CEF levels. However, since neutron pow- 
der diffraction data [7| showed that the Ce local envi- 
ronment is not much affected by the lattice distortion, 
we started fitting the low temperature spectra using the 
CEF parameters for tetragonal geometry, then added the 
effect of the molecular field of the Fc spins in the presence 
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FIG. 3: (Color online) (a) Raw S{Q,lu) spectra of CeFeAsO 
at 60 K and Ei = 30 meV on MERLIN, (b) Ce CEF ex- 
citations after subtration of the LaFeAsO background, (c) 
Temperature dependence of the 18.7 meV excitations. The 
peaks around 18 meV at 7 K are fitted with 2 lorentzians, the 
widths of which were fixed for fittings at all higher tempera- 
tures. Although the peak separation remain finite above Tm, 
the net separation caused by molecular field can be monitored 
by the change of this separation as shown in (d). A drastic 
increase of peak separation happens below 140 K, the Tm of 
Fe as shown in (e). 



of the orthorhombic structural distortion. The solid line 
in Fig. 2(a) shows our fit to the data, and the CEF pa- 
rameters for temperatures above and below Tn are also 
given in Table I. Assuming that the |I), |II), and [III) dou- 
blets in the paramagnetic state of CeFeAsO split into |1) 
to |6) singlets in the AF state as shown in Fig. 1(b), the 
observed excitations near 19 mcV should be composed of 
4 possible transitions |1) -^ |3), |1) -^ |4), |2) -^ |3), and 
|2) -^ |4). Since the transition probabilities |(1| J^z |3)| 
and |(2|J2|4)| are rather small, the observed excita- 
tions at 16.8 and 20.4 meV in Fig. 2(a) actually arise 
from |1) — *■ 1 4) and |2) — *■ |3), and are controlled by 
the thermal population of |1) and |2), respectively. This 
is consistent with the temperature dependence of these 
two excitations, where the 16.8 mode decreases and the 
20.4 meV excitation increases with decreasing tempera- 
ture [Fig. 3(c)]. Therefore, the CEF energy levels for 
|3) and |4) should be at 16.6 ± 0.7 = 17.5 meV and 20.4 
meV, respectively. 

In rare-earth substituted high- Tc copper oxides such 

as Tmo.iYo.9Ba2Cu306+a: [H and IIoo.iYo.9Ba2Cu307 

14| . the intrinsic linewidth F and peak position of the 
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FIG. 4: (Color online) Temperature dependence of the fiiox — 
18.7 meV Ce CEF excitation for CeFeAsO0.84F0.i6. (a) MER- 
LIN measurements using Ei — 30 meV at different temper- 
atures. The instrumental energy resolution is 2.1 meV at 
elastic position (horizontal bar), (b) The peak position as a 
function of temperature for the 18.7 meV CEF level. The solid 
squres are MERLIN data while open squares are FANS data, 
both show clear anomaly at Tc. (c) The intrinsic linewidth 
r(r) as a function of temperature. The solid line shows the 
expected r„(T) assuming noninateracting Fermi liquid. r(T) 
deviates from r„(r) near Tc. (d) r(r)/r„(r) shows a clear 
anomaly near Tc, consistent with (c). 



TABLE II: Wave functions of different CEF levels for 
CeFeAsO0.84F0.i6 and CeFeAsO above Tjv. 



CcFeAsOo.84Fo.: 



CeFeAsO {> Tn) 



|I) lTl/2) |Tl/2> 

|II) -0.581 |=f|) +0.814 |±|) -0.445 |t|) +0.895 |±|) 

|III) 0.814 |=Ff) +0.581 l + l) 0.895 It!) +0.446 l+l) 



rare-earth CEF excitations are used to probe the local 
magnetic response of the Cu02 planes and the formation 
of a superconducting energy gap. The CEF excitations 
of Ho and Er have also been used to study pseudogap 
and order-parameter symmetry in the underdoped super- 
conducting HoBa2Cu408 and Er2Ba4Cu70i4.92, respec- 
tively [15|. To see if the hui = 18.7 meV CEF excitation 
in CeFeAsO0.84F0.i6 is also sensitive to the occurrence 
of superconductivity, we carefully probed the tempera- 
ture dependence of its intrinsic linewidth F [Figs. 4(a), 
4(c)] and peak positions [Fig. 4(b)]. Figure 4(a) shows 
the phonon subtracted data at different temperatures ob- 
tained on MERLIN with E, = 30 meV. Figure 4(b) plots 
the peak position as a function of temperature, which 
shows a clear anomaly around Tc. Inspection of Fig. 
4(c) reveals that the linewidth decreases with decreas- 
ing temperature. For the case when the CEF levels do 
not overlap, the temperature dependence of the linewidth 



r(T) of the transition between the ground state (|1)) and 
the first excited state (|2)) in Fig. 1(c) can be accu- 
rately estimated using Eq.(l) in [ij]. To first order, this 
can be written as T{T) oc x" {Tt^i ) coth( ^^^ j ) , where 
x"(?iwi) is the local integrated magnetic susceptibility at 
fiuji (= 18.7 meV) and ks is the Boltzmann constant [13j . 
To understand r(T), we calculate the linewidth r„(T) 
below 150 K using this equation assuming a noninteract- 
ing Fermi liquid normal state [lj| , even though this is still 
under discussion. The solid line in Fig. 4(c) shows out- 
come of our calculation. The observed linewidth starts to 
deviate from the expected values near T^ as shown more 
clearly in r(r)/r„(T) [Figure 4(d)]. The most natural 
interpretation of Figs. 3(c) and (d) is that the local inte- 
grated magnetic susceptibility x" [^^i) from the Fe spin 
fluctuations at 18.7 meV increases dramatically below Tc- 
Since the temperature dependence of this excitation and 
its energy are rather similar to that of the neutron spin 
resonance recently observed in Bao.6Ko.4Fe2As2 [l8|, we 
speculate that the observed linewidth change might be 
related to the spin resonance in CeFeAsO0.84F0.i6- 

In summary, we have determined the Ce CEF levies 
in nonsuperconducting CeFeAsO and superconducting 
CcFcAsO0.84F0.i6 as a function of temperature. Our re- 
sults show that Ce CEF excitations are very sensitive to 
the Fe magnetic order in the undoped system, and dy- 
namic spin susceptibility in superconductor, and there- 
fore can be used as a probe to study spin dynamical prop- 
erties of the rare-earth based oxypnictides. 
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